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Abstract 

After more than a decade of producing a nominal 10 km aerosol product based on the 
dark target method, the MODIS aerosol team will be releasing a nominal 3 km product 
as part of their Collection 6 release. The new product differs from the original 10 km 
5 product only in the manner in which reflectance pixels are ingested, organized and 
selected by the aerosol algorithm. Overall, the 3 km product closely mirrors the 10 km 
product. However, the finer resolution product is able to retrieve over ocean closer 
to islands and coastlines, and is better able to resolve fine aerosol features such as 
smoke plumes over both ocean and land. In some situations, it provides retrievals over 
io entire regions that the 10 km product barely samples. In situations traditionally difficult 
for the dark target algorithm, such as over bright or urban surfaces the 3 km product 
introduces isolated spikes of artificially high aerosol optical depth (AOD) that the 10 
km algorithm avoids. Over land, globally, the 3 km product appears to be 0.01 to 0.02 
higher than the 10 km product, while over ocean, the 3 km algorithm is retrieving a 
is proportionally greater number of very low aerosol loading situations. Based on colloca- 
tions with ground-based observations for only six months, expected errors associated 
with the 3 km land product are determined to be greater than for the 1 0 km product: 
±0.05 ±0.25 AOD. Over ocean, the suggestion is for expected errors to be the same 
as the 10 km product: ±0.03 ±0.05 AOD. The advantage of the product is on the lo- 
20 cal scale, which will require continued evaluation not addressed here. Nevertheless, 
the new 3 km product is expected to provide important information complementary to 
existing satellite-derived products and become an important tool for the aerosol com- 
munity. 


1 Introduction 

25 The MODerate resolution Imaging Spectroradiometer (MODIS) has been observing 
the Earth from the Terra satellite since 2000 and from the Aqua satellite since 2002. 
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Among the many physical parameters derived from MODIS spectral radiances are a 
suite of products characterizing aerosol particles. There are several algorithms pro- 
ducing aerosol products from MODIS. In this paper we address the pair of algorithms 
referred to as the Dark Target algorithms, over land and ocean (Remer et al., 2005; 
5 Levy et al., 2007a, b, 2010). The original motivation behind the development of the 
Dark Target algorithms was to provide the necessary information to quantify aerosol 
effect on climate and climate processes, and thereby narrow uncertainties in estimating 
climate change (Kaufman et al., 1997, 2002; Tanre et al., 1997). Indeed in the subse- 
quent dozen years since Terra launch, the scientific literature abounds in references 
io to MODIS aerosol products to estimate direct aerosol effects (Remer et al., 2006; Yu 
et al., 2006) including the anthropogenic portion (Kaufman et al., 2005a; Christopher 
et al., 2006), to constrain climate models in their efforts to simulate climate processes 
(Stier et al., 2005; Kinne et al., 2006), to estimate intercontinental transport of aerosol 
(Kaufman et al., 2005b; Yu et al., 2012) and to further our understanding of aerosol- 
15 cloud-precipitation processes (Kaufman et al., 2005c; Koren et al., 2005; Koren and 
Wang, 2008; Loeb and Schuster, 2008). 

The fundamental scale of the MODIS Dark Target aerosol product is 10 km at nadir 
that expands roughly four-fold towards the edges of the swath. This aerosol product 
is labeled Level 2. The Level 2 data follow the orbital path of the sensor and are not 
20 gridded. Instead each retrieval is labeled by the latitude-longitude of its center. 

The Level 2 product is widely used to characterize local events, collocate with cor- 
relative data on a local level (Redemann et al., 2009; Russell et al., 2007), and to exert 
control on how the data are aggregated up to a coarser grid (Zhang et al., 2008.) Level 
2 data are automatically aggregated to a 1° x 1° global grid, labeled Level 3 and made 
25 available to the community for global-scale applications. 

One unexpected application of the MODIS Dark Target aerosol product is its use as 
a proxy for particulate pollution by the air quality community (Chu et al., 2003; Wang 
and Christopher 2003; Engel-Cox et al., 2004). The interest in using MODIS aerosol 
products to characterize air pollution has progressed both in the research arena (van 
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Donkelaar et al., 2006, 2010), and on the operational side (Al Saadi et al., 2005). The 
air quality community, except those interested in long-range and intercontinental trans- 
port of pollution (Yu et al., 2008), almost exclusively uses the Level 2 data product. Their 
interest requires identifying local areas of exposure and resolving gradients across an 
5 urban landscape. The 1 -degree global grid is insufficient. Even though the air qual- 
ity community uses the 10 km Level 2 product, there has been strong advocacy from 
that community and from others for a finer resolution product. Studies and applications 
besides air quality that would benefit from a finer resolution product include those char- 
acterizing smoke plumes from fires, those resolving aerosol loading in complex terrain 
io and those interested in aerosol-cloud processes. 

Alternative aerosol retrieval algorithms have been applied to MODIS data that pro- 
duce a finer resolution product. Most of these have been local in scope, specifically 
tuned for the local area of interest (Li, C-C. et al., 2005; Castanho et al., 2008). A 
few have been applied to a more general global retrieval over land (Hsu et al., 2004; 
15 Lyapustin et al., 2011). These finer resolution retrievals, mostly at 1 km, show much 
promise in resolving individual smoke and pollution plumes. Because there is an identi- 
fiable need for a finer resolution aerosol product from MODIS, the MODIS aerosol team 
is introducing a 3 km product as part of their Collection 6 delivery. The 3 km product will 
be a Level 2 product, available in its own files, MOD04_3K for Terra and MYD04_3K 
20 for Aqua, and offer a subset of the original parameters over both land and ocean. The 
product is created using similar structure, inversion methods and Look Up Tables as 
the basic 10 km Dark Target products. The differences arise only in the manner pixels 
are selected and grouped for retrieval. Because the MODIS Dark Target aerosol algo- 
rithms were designed with climate applications in mind and on a 1 0 km scale, they were 
25 constructed in such a way to suppress noise in the retrieval. The danger of applying a 
similar inversion scheme on a finer scale is the possibility of introducing noise. 

In this paper we introduce the MODIS Dark Target 3 km product by reviewing the 
algorithm producing the 1 0 km product and detailing the changes that allow for retrieval 
at 3 km. Then we demonstrate the new product in side-by-side comparisons between 
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the 3 km and 10 km retrievals of the same scenes. Finally, we produce a limited vali- 
dation of the new product based on collocations with ground-based sunphotometry on 
a global basis, but for only 6 months of Aqua data. A companion paper in this same 
special issue describes the application of the new product with greater detail across an 
5 urban/suburban landscape (Munchak et al., 2012). 


2 MODIS aerosol retrieval at 10 km and 3 km resolution 

The MODIS Dark Target aerosol algorithms are well-documented in the literature (Kauf- 
man et al., 1997; Tanre et al., 1997; Remer et al., 2005, 2012; Levy et al., 2007a, b, 
2012). Here we provide a short review in order to highlight the differences between the 
io 10 km and 3 km algorithms. 

The MODIS Dark Target aerosol algorithms are two separate algorithms, one ap- 
plied over ocean and one over land. They operate on five-minute segments of MODIS 
along-orbit data, known as “granules”. Over ocean, the inputs consist of the MODIS- 
measured geolocated radiances normalized to reflectance units in 7 wavelengths (0.55, 
is 0.66, 0.86, 1.24, 1.38, 1.63 and 2.13 pm), total column ozone concentrations from the 
NOAA Office of Satellite Product Operations, total column precipitable water vapor 
from the National Center for Environmental Prediction (NCEP) reanalysis, the MODIS 
cloud mask (MOD/MYD35) and in Collection 6 the surface wind speed from NCEP. The 
over land algorithm uses the 0.47, 0.66, 0.86, 1.24, 1.38, 2.13 pm channels, the ozone 
20 concentrations, and the total column precipitable water vapor. The 1 .38 pm channel is 
available at 1 km resolution and is used to mask clouds, not retrieve aerosol. The 0.66 
and 0.86 pm channels are available at 0.25 km resolution, and the other channels are 
available at 0.5 km resolution. Over land the 0.66 and 0.86 channels are used at their 
native resolution to identify inland and ephemeral water sources, then these bands 
25 are averaged to 0.5 km to create collocated spectral data at 0.5 km resolution, over 
both land and ocean. Normally, one granule is composed of 2708 x 4060 pixels at this 
resolution. 
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2.1 The MODIS 10 km Dark Target aerosol retrieval 

Before the MODIS pixels are organized into retrieval boxes gaseous correction is ap- 
plied and individual pixels are marked for either the ocean retrieval or the land retrieval. 
Then, for the 10 km (nominal at nadir) retrieval, we organize the entire MODIS granule 
5 into groups of 20 x 20 pixels, which we refer to as “retrieval boxes”. The left side of 
Fig. 1 illustrates a 10 km retrieval box outlined in magenta. 

After organization, the first step is to select the pixels within the retrieval box to be 
used in the retrieval. The algorithm avoids clouds, ocean sediments, glint, snow, ice, in- 
land water and bright surfaces. Clouds are identified by means of spatial variability, ratio 
io and threshold tests with additional assistance from specific tests from the MOD/MYD35 
product (Martins et al., 2002; Gao et al., 2003; Frey et al., 2008; Remer et al., 2012). 
Sediments are identified in the ocean using spectral tests (Li et al., 2003) and sun glint 
is eliminated through the use of a 40° glint mask. Snow and ice are identified using 
spectral tests (Li, R-R. et al., 2005). Subpixel inland water is identified using the 0.66 
is and 0.86 pm channels and bright surfaces (> 0.25) at 2.13 pm are avoided altogether. 
After all unsuitable pixels are identified and deselected, the remaining pixels are sorted 
from lowest to highest reflectance at 0.86 pm over ocean and at 0.66 pm over land. The 
darkest and brightest 25% of remaining pixels in the retrieval box are arbitrarily des- 
elected over ocean, and the darkest 20% and brightest 50% of the remaining pixels 
20 are deselected over land. This means that in a 20 x 20 box, there are at most, 200 pix- 
els over ocean or 120 pixels over land. In a 10 km retrieval box with 400 pixels even if 
many pixels are avoided for one or more of the above reasons or arbitrarily deselected 
at the dark or bright end of the distribution, there may still exist sufficient uncontami- 
nated pixels to represent aerosol conditions in that box (Remer et al., 2012, this issue). 
25 The ocean algorithm requires 10 out of the 400 pixels at the 0.86 pm channel and at 
least 30 pixels total distributed across channels 0.55, 0.66, 1.24, 1.63 and 2.13 pm to 
represent aerosol conditions and to produce a high quality retrieval in that box. The 
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land algorithm requires 51 pixels in the 0.66 pm channel, for a high quality retrieval, but 
only 12 pixels for a degraded quality retrieval. 

After the selection procedure if there are sufficient pixels remaining, the mean spec- 
tral reflectance is calculated from the remaining pixels. From the 400 pixels in the re- 
trieval box, there emerges a single set of spectral reflectance representing the cloud- 
free conditions in the box. The spectral reflectances are matched in a Look Up Table 
with pre-calculated values. Over ocean the entries in the Look Up Table are calcu- 
lated assuming a rough ocean surface, and in Collection 6 there are separate Look Up 
Tables for different surface wind speeds as determined from the NCEP wind data (Klei- 
dman et al., 2012; Levy et al., 2012). Over land, the surface reflectance is constrained 
by spectral functions relating the visible to 2.13 pm (Levy et al., 2007a). 

TJ 

2.2 The 3 km retrieval 

0 

“"S 

The only differences between the 3 km algorithm and the 1 0 km algorithm are the way 
the pixels are organized and the number of pixels required to proceed with a retrieval 
after all masking and deselection are accomplished. Figure 2 presents a flow chart 
showing the separate paths for the 10 km and 3 km retrievals. The black boxes run- 
ning along the center of chart identify processes that are identical in both retrievals. 

The inputs are identical, as are the masking procedures. The exact same 0.5 km pix- 
els identified as cloud, sediment etc. in the 10 km algorithm are identified as cloud, 
sediment etc. in the 3 km algorithm. The difference is in how the two algorithms make 
use of these 0.5 km designations. Once the 3 km algorithm has identified the pixels 
suitable for retrieval and decided that a sufficient number of these pixels remain, the 
spectral reflectances are averaged and the inversion continues exactly the same as in 
the 10 km algorithm. The same assumptions are used, the same Look Up Tables, the 
same numerical inversion and the same criteria to determine a good fit. 

In the 3 km retrieval the 0.5 km pixels are arranged in retrieval boxes of 6 x 6 arrays 
of 36 pixels, illustrated by the schematic in the right hand side of Fig. 1 . Note that in 
the 3 km retrieval box, the exact same pixels identified as cloudy in the 10 km retrieval 
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box (denoted by the white rectangles) are identified as cloudy in the 3 km box. This 
is because both algorithms apply identical criteria to masking undesirable pixels. The 
3 km retrieval attempts to apply similar deselection of pixels at the darkest and bright- 
est ends of the distribution: 25% and 25% over ocean, and 20% and 50% over land. 
5 Once these darkest and brightest pixels are discarded, the algorithm averages the re- 
maining pixels to represent conditions in the 3 km retrieval box. The algorithm requires 
a minimum of 5 pixels at 0.86 pm over ocean with at least 12 pixels distributed over the 
other five channels and 5 pixels are required over land in order to continue and make a 
retrieval. This is actually a more stringent requirement for ocean (1 4 % of 36), than what 
io is required by the 1 0 km retrieval (2.5 %) for the best quality retrieval. The requirement 
over land is about the same in the 3 km retrieval as it is in the 1 0 km retrieval (14 % and 
13%, respectively). 

Tables 1 through 3 list the parameters available in the MOD04_3K and MYD04_3K 
files. Not all of the diagnostics available at 10 km are included at 3 km, but 
is most of the parameters are there. The data set includes an integer Quality Flag 
(Land_Ocean_Quality_Flag) that designates each 3 km retrieval as “3”, “ 2 ”, “ 1 ” or “ 0 ”. 
The same recommendations apply to the 3 km product as to the original product. 
Ocean retrievals are valid for all nonzero Quality Flags while land retrieval products are 
only recommended for Quality Flags = “3”. The 3 km product also includes more de- 
20 tailed diagnostics about the retrieval embedded in the Quality .Assurance .Ocean and 
Quality _Assurance_Land parameters. Note, that the criteria to fill these quality diagnos- 
tics differ slightly from the criteria used to fill the same-named parameters at 10 km. 
Detailed information on the quality diagnostics will be available in the Collection 6 
Algorithm Theoretical Basis Document available at http://modis-atmos.gsfc.nasa.gov/ 
25 reference_atbd.html. 

The 3 km product is designed to provide insight into the aerosol situation on a fo- 
cused local basis and will not be aggregated to the Level 3 global 1 -degree grid. All 
Level 3 MODIS aerosol products will be derived from the 10 km product. 
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The 3 km retrieval should closely mirror the results of the 1 0 km retrieval because the 
majority of the two algorithms are identical, but it should be able to resolve gradients 
across the 10 km retrieval box that would otherwise be missed. There is a possibility 
that the new algorithm will introduce additional noise, especially over land where pixels 
5 representing inhomogeneous surfaces would have been eliminated in the deselection 
process of the larger box, but are now included in the retrieval. However, this tendency 
towards noise may be mitigated by the slightly more stringent requirements in deselec- 
tion and the minimum number of pixels needed to represent the box. The advantageous 
mitigation features will be more apparent over ocean than over land. Overall, the two 
io algorithms provide different sampling of the aerosols over the globe and because of 
this alone, are not expected to provide the same global statistics. 


3 Examples of results from 3 km algorithm 

The 3 km algorithm was applied to six months of MODIS data from Aqua-MODIS: Jan- 
uary and July 2003, 2008 and 2010. This is a special database of MODIS data used 
15 to test new MODIS algorithms before implementation into operational production. The 
input radiances and MOD/MYD35 cloud mask are not final versions that appear or 
will appear in a publicly available data Collection, but the inputs are traceable within 
the MODIS processing environment. In Figs. 3-6, we show examples from Day 183, 
15 July 2008 and from Day 12, 12 January 2010 that illustrate the new 3 km product 
20 and how it differs from the 10 km product applied to exactly the same input data. 

Figure 3 compares the 10 km and 3 km ocean retrievals at over the Mediterranean 
Sea off the coast of Tunisia and Libya during a moderate dust event. The two resolu- 
tion products produce almost the exact same aerosol field with the same gradient and 
same magnitude aerosol optical depth. This is because the two algorithms are essen- 
25 tially the same. The only difference is that the finer resolution product is able to make 
retrievals closer to the small islands in the image. We find that this is typical of the 
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3 km product. It offers over-ocean retrievals closer to land, nearer to islands and within 
narrow waterways and estuaries. 

Figures 4 illustrate the apparent advantage of the 3 km product to resolve smoke 
plumes from fires. The fire is a large wild fire burning in Canada. The 10 km product 
5 does not capture the long narrow smoke plume leading towards the northwest, but 
the 3 km product does. One of the major advantages of the 3 km product is its ability 
to better resolve smoke plumes than the 10 km product. Even so, because the cloud 
identification algorithm in the 3 km product is the same as in the 10 km product, based 
primarily on spatial variability, the 3 km product still improperly confuses the thickest 
io parts of the smoke plume with a cloud and mistakenly refuses to retrieve there. 

Figure 5 demonstrates the potential for different sampling by the two products. The 
situation is a highly polluted episode over much of southeastern China. Here the 3 km 
algorithm makes retrievals over a broad area, while the 10 km algorithm finds few op- 
portunities to retrieve. The few places of overlap result in similar values of aerosol 
15 optical depth. The only AERONET station in the image is at Hong_Kong_PolyU (22°1 8', 
1 14°1 T), which reports a collocated AOD interpolated to 0.55 pm at MODIS overpass 
time of 0.38. The 10 km algorithm does not produce a retrieval at this station, but the 
3 km algorithm does, producing an AOD of 0.45, a reasonable match. The collocation 
procedure and quantification of expected uncertainties are described in Sect. 5 below. 
20 Figure 6 shows a potential concern of switching indiscriminately to the 3 km product. 
In this retrieval over the highly urbanized surface of Los Angeles and environs, the sur- 
face is incompatible with the current version of the Dark Target retrieval. The seasoned 
pixel selection process of the 10 km algorithm is able to recognize this incompatibility 
and chooses not to retrieve over Los Angeles. However, the 3 km product does retrieve, 
25 and the result is a scattering of retrieved AOD > 0.8 over the region. Although there is 
no ground truth to determine whether these points are accurate high AOD situations or 
noisy artifacts of the retrieval, it is highly likely that they are artifacts that the 3 km re- 
trieval fails to avoid. Although the results of the 3 km product mirror the 1 0 km retrievals, 
we do find an increase of noisy artifacts in the finer resolution product. Unfortunately 


o 


( / ) 
O 
c 

V) 

(f) 

o' 


AMTD 

6, 69-112, 2013 


TJ 


CD 

"O 

CD 


D 

c / )' 
o 
c 

C/) 

c n 
o' 


Tl 


CD 

"O 

CD 


MODIS 3 km aerosol 
product: algorithm 
and global 
perspective 

L. A. Remer et al. 


Title Page 


D 

< ri 
o 
c 

C/) 

cn 

o' 

o 

"0 

CD 

■o 

CD 


D 

c/)' 

o 

c 

( f) 

c n 
o' 

D 

"0 

CD 

“O 

CD 


Abstract 


Introduction 



Conclusions 


References 



Tables 


Figures 



Printer-friendly Version 
Interactive Discussion 



78 

















this occurs most frequently over urban surfaces, a type of location of most interest to 
the air quality community. 

c n 
o' 

D 

4 Global mean aerosol statistics of the 3 km product 

■ o 
CD 

The global mean aerosol optical depth calculated from the 10 km and 3 km algorithms 
are different because the two algorithms sample differently. In general the global statis- 
tics of the 3 km product tracks the day-to-day variation of the 10 km product. Figure 7 
shows the day-to-day global mean AOD differences between the two products over 
ocean and land for the three months of January merged into one continuous time se- 
ries for plotting purposes and the three July months merged into another. A positive 
difference indicates that the 3 km AOD is higher than the 10 km AOD. The differences 
in global monthly mean AOD between the two resolution data sets is -0.004 in January 
and nearly 0.000 in July over ocean and 0.003 in January and 0.010 in July over land. 

The largest day to day differences between the products is seen over land in January, 
where daily differences can be either positive or negative. Because of snow cover over 
the northern land masses in January there are relatively fewer retrievals contributing 
to the global mean AOD, causing relatively larger day-to-day fluctuations over land in 
this month as compared with July. The large daily difference in the products in January 
over land, corresponding to Day 12 on 2010 on the graph (12 January 2010) does not 
correspond to an unusual spike in global land AOD, just to a large difference between 
the two products. On this day the differences are concentrated in a handful of granules, 
of which the most dramatic is the 05:35 UTC granule shown in Fig. 5. Likewise other 
spikes on the difference plots are not associated with unusual global mean AOD. 

Figure 8 shows the scatter plots of the 3 km daily global mean over land and ocean 
plotted against the same using the 10 km product. Januarys and Julys are denoted by 
different symbols. The correlation between the two resolutions over land in July is very 
high ( R = 0.997) with a regression equation slope approaching 1 and offset of only 
0.015. There is greater deviation between the two resolutions in global mean statistics 
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on a day-by-day basis over land in January and over ocean in both months, but the 
correlation remains high ( R > 0.95). January land is noisiest with the lowest correlation, 
most likely due to the snow cover and senescent vegetation of the northern temperate 
latitudes that severely reduces the number of retrievals on any day. As we see in Fig. 5, 
5 the two resolution algorithms are affected differently by these less than optimal surface 
conditions resulting in significantly different samples from which to construct the global 
mean. The result is the lower correlation of Fig. 8. Over ocean in January, the 3 km 
product tends to produce lower global mean AOD than the 10 km product. In July the 
tendency for lower global mean AOD is seen in the lower AOD range, but there is a 
io strong positive slope in the regression equation so that when 10 km AOD is greater 
than 0.12, the 3 km product actually tends towards higher AOD. As of now we do not 
have an explanation for these tendencies in the global mean AOD statistics over ocean. 

Figure 9 shows the histograms calculated from the six months of data described in 
Sect. 3, with all Januarys combined and all Julys combined, land and ocean separately, 
is These histograms are constructed from individual retrieval boxes accumulated for the 
entire three-month periods with no spatial or diurnal averaging. The histograms are 
plotted with relative frequency rather than total number of retrievals in each bin because 
the 3 km product at finer spatial resolution produces approximately 1 1 times and 7 
times the number of retrievals produced by the 10 km product, over land and ocean, 
20 respectively. Here we see over land with the 3 km product a decrease in the proportion 
of negative and very low AOD retrievals and an increase in retrievals in the 0.05- 
0.15 range in January and in the 0.15-0.35 range in July. This shift is more apparent in 
Fig. 10. The overall increase in global mean AOD over land with the 3 km product noted 
in Figs. 7 and 8 appears to be due to shifting to moderate AOD from very low AOD and 
25 not from introducing spikes at very high AOD. In contrast, over ocean the histograms 
show an increase in very low AOD (0-0.05) at the expense of slightly higher AOD 
(0.05-0.15 or 0.20). Figure 10 clarifies the differences between the 3 km and 10 km 
histograms. 


o 


( / ) 
O 
c 

V) 

(f) 

o' 


AMTD 

6, 69-112, 2013 


TJ 


CD 

"O 

CD 


D 

c / )' 
o 
c 

C/) 

c n 
o' 


Tl 


CD 

"O 

CD 


MODIS 3 km aerosol 
product: algorithm 
and global 
perspective 

L. A. Remer et al. 


Title Page 


D 

< ri 
o 
c 

C/) 

cn 

o' 

o 

"0 

CD 

■o 

CD 


D 

c/)' 

o 

c 

( f) 

c n 
o' 

D 

"0 

CD 

“O 

CD 


Abstract 


Introduction 



Conclusions 


References 



Tables 


Figures 



Printer-friendly Version 
Interactive Discussion 



80 

















Overall the 3 km product tracks the 1 0 km product on a day-to-day basis, although 
the product over land tends to be higher than the 10 km product, and over the ocean 
lower. There is seasonal variation to these tendencies. 

D 

TJ 

Q) 

5 Global validation 

— s 

The six months of test data described above (January and July 2003, 2008 and 2010) 
were collocated with Level 2.0 AERONET observations to test the accuracy of the 
retrievals. For the 10 km retrievals we use the Petrenko et al. (2012) protocol for collo- 
cations, which differs slightly from the one introduced in Ichoku et al. (2002). Here, a 
collocation is the spatio-temporal average of all AERONET AOD measurements within 
±30 min of MODIS overpass and the spatial average of all MODIS retrievals within a 
25 km radius around the AERONET station. At nadir, the 25 km radius can encompass 
roughly 25 MODIS retrieval boxes, each at 10 km. However, MODIS pixel resolution 
increases with scan angle, as does the size of the retrieval boxes. At swath edges, 
the aerosol product box can be approximately 40 km instead of 10 km. The collocation 
protocol still calls for a 25 km radius, which now encompasses only parts of 4 boxes. 

To be included in the analysis, a 10 km collocation must include two AERONET obser- 
vations within the hour and at least 20% of the potential MODIS retrievals. This would 
be at least 5 out of 25 possible retrievals at nadir, but only 1 out of a possible 5 re- 
trieval boxes towards the scan edge. For the 3 km retrievals we apply a 7.5 km radius 
around the AERONET stations, which encompasses 25 MODIS 3-km retrieval boxes 
at nadir. Again, for 3 km retrievals, we require at least 2 AERONET observations and 
20% of possible MODIS retrievals for the collocation to be included in the analysis. 

Note that the collocations are filtered using MODIS Quality Assurance (QA) flags. Only 
those MODIS retrievals with QA = 3 over land and QA > 0 over ocean within the 25 km 
or 7.5 km radius are included in the statistics of the collocation and the requirement 
is for at least 20% of MODIS retrievals with acceptable QA for the collocation to be 
included in the analysis. The only wavelength examined is at 550 nm. This requires the 
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AERONET values to be interpolated to this wavelength in order to match MODIS. A 
quadratic fit in log-log space is used to make the interpolation (Eck et al., 1999). 

Figure 1 1 shows the binned scatter plots from six months of collocations for land. 
There are 3298 collocations at 10 km and 3283 at 3 km. The data were sorted ac- 

TJ 

cording to AERONET AOD and bins designated for every 50 collocations. The mean 
and standard deviation of both the AERONET and MODIS AODs were calculated. The 
mean values for each bin are plotted and the error bars indicate ±1 standard deviation. 

The red line represents the linear regression calculated for the full data set of over 
3000 collocations and not the 66 bins plotted in the figure. The dashed lines indicate 
the “expected error” (Levy et al., 2010). For land the expected error is 

c n 

AAOD = ±0.05 ±0.15 AOD 

TJ 

0 ) 

Table 4 provides the mean AOD of each data set, the correlation coefficient, the regres- 
sion statistics and the number of MODIS retrievals that fall within expected error, fall 
above the expected error bound (the upper dashed line) and fall below the expected 
error bound (the lower dashed line in Fig. 11). 

Figure 12 shows similar binned scatter plots for the ocean retrieval. There are fewer 
ocean collocations, 1100 at 10 km and 697 at 3 km. Because AERONET stations are 
generally on land, the 10 km retrieval with its longer radius will intercept more ocean 
retrievals than the 3 km retrieval with only a 7.5 km radius. The expected error for ocean 
is 

AAOD = ±0.03 ± 0.05 AOD. 

D 

These figures show that the 10 km product of this six-month database is meeting its 
expectations with most retrievals falling within expected error (71 % over land and 69% 
over ocean). Correlations are high (0.87 over land and 0.93 over ocean) and most re- 
trievals and the linear regression fall close to the 1:1 line. The land retrieval does show 
a positive bias of ~ 0.015 and the ocean ~ 0.005, but these are well-within expected 
uncertainties. 
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The 3 km product is also highly correlated to AERONET observations and most of 
the 3 km retrievals fall within the expected error bounds. However, the 3 km product 
matches AERONET less well than does the 10 km MODIS product. The most obvious 
degradation of accuracy between the 3 km and 10 km products is over land where the 
5 finer resolution product has more than doubled its positive bias. This can be seen in 
Fig. 1 1 where the red cluster of points and regression line are now above the green 
1:1 line. The ocean retrieval has also significantly increased its positive bias against 
AERONET at finer resolution. In addition to the increased bias the 3 km product intro- 
duces additional noise at low AOD. This can be seen in Fig. 12, where bin standard 
io deviations have increased, especially in the two lowest bins. 

The positive bias in the MODIS-AERONET collocations is best seen in Fig. 13, where 
the difference in each bin has been plotted against AERONET AOD, separately for 
each product resolution, land and ocean. The x-axis in the figure has been truncated 
at AERONET AOD = 0.6 to emphasize the lower AOD range. Over ocean, the lowest 
is bins show a large relative increase in AOD bias, but for most of the AOD range the 
3 km and 1 0 km biases against AERONET are approximately the same. Over land, the 
positive bias in the 3 km AOD exceeds 0.05 from AERONET AOD = 0.15 to 0.50. 

Because AERONET stations used in the ocean validation collocations are limited 
to coastal and island stations they do not provide adequate sampling over the open 
20 ocean where most of the ocean retrievals take place. The Maritime Aerosol Network is 
a relatively new network under the AERONET umbrella that offers quality-assured AOD 
measurements from ships of opportunity over the open ocean (Smirnov et al., 2009). 
We identified collocations between the MAN AOD and the Aqua-MODIS 10 km and 
3 km retrievals during the six months of this test data set, January and July, 2003, 2008 
25 and 2010. A collocation was defined using the same spatial averaging and minimum 
number of MODIS retrieval criteria as with a MODIS-AERONET collocation. The MAN 
measurements temporally average a series of measurements together to form one 
measurement. All “series” measurements within a half hour of the MODIS overpass are 
temporally averaged. Due to the limited amount of MAN data, there was no minimum 


o 


( / ) 
O 
c 

V) 

(f) 

o' 


AMTD 

6, 69-112, 2013 


TJ 


CD 

"O 

CD 


D 

c / )' 
o 
c 

C/) 

c n 
o' 


Tl 


CD 

"O 

CD 


MODIS 3 km aerosol 
product: algorithm 
and global 
perspective 

L. A. Remer et al. 


Title Page 


D 

< ri 
o 
c 

C/) 

cn 

o' 

o 

"0 

CD 

■o 

CD 


D 

c/)' 

o 

c 

( f) 

c n 
o' 

D 

"0 

CD 

“O 

CD 


Abstract 


Introduction 



Conclusions 


References 



Tables 


Figures 



Printer-friendly Version 
Interactive Discussion 



83 

















number of MAN “series” measurements required to create a collocation. Only ocean 
retrievals with a QA > 1 are included in the collocation. Because of the requirement that 
there be at least 20 % of available retrievals within a 7.5 km radius for the 3 km product 
and a 25 km radius for the 1 0 km product, the two sets of collocations do not necessarily 
5 contain the same points. In our sets of MAN-MODIS collocations both the 3 km and 
10 km retrievals contain 37 points, of which 27 are the same and 10 are unique for 
each resolution. The AERONET collocation exercise also results in differently sampled 
data sets, but the results of this different sampling is much more apparent in the MAN 
process because of the small statistical sampling size. 
io The results of the MODIS-MAN collocations are shown in Fig. 14. Most MAN collo- 
cations occurred in very clean environments with the exception of 3 points (only 2 at 
10 km) that occurred in the tropical Atlantic in the dust belt and a moderate value near 
AOD ~ 0.2 from the Hudson Bay that appears to be wildfire smoke. The highest AOD 
collocation in the 3 km data base (9 July 2008 from the dust belt of the Atlantic Ocean 
15 with MODIS AOD = 0.67) does not appear in the 1 0 km data base. Most of the MODIS 
3 km points agree with the MAN observations to within expectations. The exceptions in- 
clude that one high dust point and a group of collocations at low AOD that are all traced 
to a specific cruise within sight of the Antarctic coast in January 2008. The regression 
and correlation statistics are given in Fig. 14 and in Table 4 but because of the outlying 
20 dust point in the 3 km set, these statistics are not robust descriptions of the bulk of the 
observations. For example while the mean AODs of the MAN and MODIS populations 
are 0.077 and 0.101, respectively, the medians are 0.051 and 0.059, respectively. The 
37 points are inadequate to fully represent the relationship between MODIS 3 km re- 
trievals over the wide variety of conditions experienced over the world’s oceans. They 
25 are shown here to supplement the other inadequate data set of the AERONET coastal 
and island sites. Together the available ocean validation does suggest, without firm 
proof, that the ocean 3 km retrieval will achieve similar levels of uncertainty as the well- 
studied ocean 10 km product. 
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6 Discussion and conclusions 

The MODIS Dark Target aerosol algorithm relies on a data selection process that iden- 
tifies a relatively few ideal pixels to use in the retrieval of aerosol optical depth and 
other aerosol characteristics. By choosing only a few pixels to represent aerosol over 
5 a moderate resolution retrieval box, noise is reduced and situations difficult to retrieve 
are avoided. Inherent in this selection procedure is an assumption that aerosol proper- 
ties do not vary across the retrieval box, so that the aerosol conditions across the box 
can be represented by just a small fraction of pixels. Except for specific situations near 
sources: smoke plumes from fires, dust plumes from playas, etc., aerosol homogeneity 
10 over mesoscale lengths of 40-400 km has been considered to be a robust assumption 
(Anderson et al., 2003). The 10 km (nadir) to 40 km (swath edge) retrieval is a reason- 
able algorithm construct, given this understanding of aerosol homogeneity. However, as 
our opportunities to observe aerosols increase and our understanding grows, we know 
now that aerosol may vary frequently over much smaller spatial scales (Shinozuka and 
is Redemann, 201 1 ; Munchak et al., 201 2 this issue). Not only will the 1 0 km retrieval box 
lose the details of local variability, the assumption on which it is based may be in error. 

Because there is need for finer resolution aerosol products to resolve individual 
plumes and fine gradients, the MODIS Science Team is introducing a 3 km product in 
their Collection 6 delivery. This product will be available at the granule level, in separate 
20 files labeled as MOD04_3K and MYD04_3K, for Terra and Aqua, respectively. The new 
product differs from the original 10 km product only in the manner in which reflectance 
pixels are ingested, organized and selected by the aerosol algorithm. All cloud, sur- 
face, sediment, snow and ice masking remain identical to the original algorithm, and 
Look Up Tables and inversion methods have not changed. The only difference is in 
25 how the algorithm arbitrarily discards additional good pixels to obtain the best pixels 
for retrieval. In the 3 km algorithm, pixels will be used for retrieval that would have been 
arbitrarily discarded by the 10 km algorithm. 
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The 3 km product exhibits expected characteristics. It resolves aerosol plumes and 
details of fine-scale gradients that the 10 km product misses. In some situations it pro- 
vides retrievals over entire regions that the 10 km barely samples. The 3 km product 
also allows the ocean retrieval to retrieve closer to islands and in narrow bays. On the 
other hand, in situations known to be difficult for the Dark Target retrieval, such as over 
bright surfaces and especially over urban surfaces, the 3 km retrieval introduces spo- 
radic unrealistic high values of AOD that are avoided more successfully by the 10 km 
retrieval. We can label these artifacts as “noise”, but it is not random noise because the 
tendency over land is to over estimate AOD in these artifacts. 

Over land, globally, the 3 km product appears to be 0.01 to 0.02 higher than the 
10 km product. There are strong differences between January and July, indicative of a 
seasonal shift, but with only 6 months of data over 3yr to analyze, the seasonal pattern 
cannot be resolved. The fact that over land AOD is higher in the 3 km product than in 
the 10 km product could be due to the fact that the finer resolution product retrieves in 
strong aerosol plumes, missed by the coarser resolution product. This would suggest 
that the 3 km global mean AOD is a truer representation of the global mean. On the 
other hand, the 3 km product introduces a high-biased noise over bright and/or urban 
surfaces, and so, the global mean from the 10 km product would remain the truer repre- 
sentation. Collocations with AERONET observations suggest the latter. The 3 km AOD 
over land compares less well with AERONET than does the 10 km AOD, decreasing 
correlation, increasing high bias and shifting retrievals from within expectations of un- 
certainty to exceeding those expectations. We conclude that the 3 km AOD over land is 
less accurate and less robust than the 10 km AOD. Estimated uncertainty of the 3 km 
land product, based on 6 months of collocations with AERONET, suggest that 67% of 
retrievals should fall within 

< f) 

AAOD_3 km = ±0.05 ± 0.25 AOD 

"O 

with the understanding that most retrievals will fall within the positive end of this error 
bounding, leaving a positive bias. 
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Over ocean, globally, the 3 km algorithm picks up proportionally a greater number 
of very low AOD cases than the 10 km algorithm. Again, just from comparing the two 
MODIS resolution products we cannot tell whether this low bias is a better representa- 
tion of global AOD or not. Unfortunately, the comparison with AERONET cannot either. 
5 AERONET stations used for the over ocean validation are isolated to a limited number 
of island and coastal locations. The very low AOD situations tend to occur over open 
ocean. The AERONET analysis suggests that the 3 km algorithm introduces positive 
bias, not negative. This is because the 3 km product retrieves closer to shore, where the 
AERONET station is located. In these locations, likelihood of sediment contamination 
io is high, aerosol is continental in nature and we expect that the result of the AERONET 
validation over ocean is not applicable to the global oceanic AOD retrieval. We can- 
not even conclude about the relative accuracy of the two products in the coastal zone, 
because the validation procedure enables the 10 km product to encompass retrievals 
much further from shore than the 3 km product. However, in the coastal zone, within 
15 7.5 km of shore, we conclude that 69 % of 3 km retrievals should fall within 

AAOD_3 km = ±0.03 ± 0.05 AOD 

with caution expressed for the lowest AOD situations, which appear to be randomly 
noisier than the expressed uncertainty bound. 

The data from the MAN cruises were also inadequate to state firm conclusions about 
20 the 3 km ocean retrieval because only 37 collocations were identified in the analysis 
data set, and several of those points were highly localized to specific locations such 
as Hudson Bay and near the Antarctic coast. However, even in this limited data set, 
68 % of the 3 km retrievals were contained within the error bounds stated above for the 
ocean retrieval. 

25 All analysis presented in this paper represent Aqua-MODIS from a limited 6 month 
analysis data set. Terra-MODIS data from the same 6 months were also examined sub- 
sequently, and no results from the Terra analysis contradict the conclusions presented 
here. 
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We have made no attempt in this analysis to address cloud effects on the 3 km 
retrieval. Some of the high bias seen in the over land product could be due to additional 
cloud contamination and cloud effects in the product instead of artifacts introduced by 
bright surfaces. However, the fact that the ocean 3 km product does not also contain this 
high bias prompts us towards considering surface effects and not clouds, but without 
further analysis we cannot make firm conclusions. Definitely future work will need to 
address cloud 3-D effects (Wen et al., 2007; Marshak et al., 2008), the so-called twilight 
zone or continuum (Koren et al., 2007; Charlson et al., 2007) and other cloud issues 
on the finer resolution aerosol product. 

Overall the 3 km product mimics the 10 km product, globally, and on a granule-by- 
granule basis. Because the new product is essentially the same as the traditional Dark 
Target product, it is well-understood and the limited analysis presented above is suffi- 
cient to recommend its use by the community with the following caveats: 

- Global studies should continue to make use of the more robust and well-studied 
10 km product. The 3 km product’s use should be restricted to obvious situations 
that require finer resolution. 

- Only the AOD at 550 nm was examined in this study. Differences in the spectral 
AOD and size parameter retrievals over ocean in the two resolution products are 
possible. 

- Aerosol-cloud studies with the 3 km product should proceed cautiously. At this 
time, we do not know specifically how the 3 km product is affected in the proximity 
of clouds. 

- While the air quality community will be eager to apply the 3 km product across an 
urban landscape, this must proceed cautiously because of known artifacts in the 
product over urban surfaces. See Fig. 6 and Munchak et al. (2012 this issue). 

The power of the new product is on the local scale, not the global one, as was stud- 
ied here. Future work that applies the MODIS 3 km aerosol product to local aerosol 
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situations in case studies and evaluates the results will be necessary to continue the 
work started here (i.e. Munchak et al., 2012). We expect the new 3 km product to pro- 
vide important information complementary to existing satellite-derived products and 
become an important tool for the aerosol community. 
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Table 1. Land and ocean parameters of the MOD/MYD04_3K file and the parameter’s dimen- 
sion. 


Parameter 

dimension 

Longitude 

(X,Y) 

Latitude 

(X,Y) 

Scan_Start_Time 

(X,Y) 

Solar_Zenith 

(X,Y) 

Solar_Azimuth 

(X,Y) 

Sensor_Zenith 

(X,Y) 

Sensor_Azimuth 

(X,Y) 

Scattering_Angle 

(X,Y) 

Glint_Angle 

(X,Y) 

Land_Ocean_Quality_Flag 

(X,Y) 

Land_sea_Flag 

(X,Y), 

0ptical_Depth_Land_And_0cean 

(X,Y) 

lmage_0ptical_Depth_Land_And_0cean 

(X,Y) 


(X,Y) refers to a 2-dimensional array along and across the swath. 
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Table 2. Ocean parameters of the MOD/MYD04_3K file and the parameter’s dimension. 


Parameter 

dimension 


Wind_speed_Ncep_Ocean 

(X,Y) 


Solution Jndex_Ocean_Small 

(X,Y) 


Solution_lndex_Ocean_Large 

(X,Y) 


Effective_0ptical_Depth_Best_0cean 

(X,Y,7); 0.47, 0.55, 0.66, 
1.24, 1.63, 2.13 pm 

CD 

00 

o 

Effective_0ptical_Depth_Average_0cean 

(X,Y,7); 0.47, 0.55, 0.66, 
1.24, 1.63, 2.13 pm 

c 6 
00 
o 

Optical_Depth_Small_Best_Ocean 

(X,Y,7); 0.47, 0.55, 0.66, 
1.24, 1.63, 2.13 pm 

co" 

00 

o 

Optical_Depth_Small_Average_Ocean 

(X,Y,7); 0.47, 0.55, 0.66, 
1.24, 1.63, 2.13 pm 

cd" 

00 

o 

Optical_Depth_Large_Best_Ocean 

(X,Y,7); 0.47, 0.55, 0.66, 
1.24, 1.63, 2.13 pm 

CD 

00 

o 

Optical_Depth_Large_Average_Ocean 

(X,Y,7); 0.47, 0.55, 0.66, 
1.24, 1.63, 2.13 pm 

c 6 
00 
o 

Mass_Concentration_Ocean 

(X,Y) 


Aerosol_Cloud_Fraction_Ocean 

(X,Y) 


Effective_Radius_Ocean 

(X,Y,2): best, average 


PSML003 

(X,Y,2): best, average 


Asymmetry_Factor_Best_Ocean 

(X,Y,7); 0.47, 0.55, 0.66, 
1.24, 1.63, 2.13 pm 

c 6 
00 
o 

Asymmetry_Factor_Average_Ocean 

(X,Y,7); 0.47, 0.55, 0.66, 
1.24, 1.63, 2.13 pm 

c 6 
00 
o 

Backscattering_Ratio_Best_Ocean 

(X,Y,7); 0.47, 0.55, 0.66, 
1.24, 1.63, 2.13 pm 

co" 

00 

o 

Backscattering_Ratio_Average_Ocean 

(X,Y,7); 0.47, 0.55, 0.66, 
1.24, 1.63, 2.13 pm 

co" 

00 

o 

Angstrom_Exponent_1 _Ocean (0.55/0.86 micron) 

(X,Y,2): best, average 


Angstrom_Exponent_2_Ocean (0.86/2.1 micron) 

(X,Y,2): best, average 


Least_Squares_Error_Ocean 

(X,Y,2): best, average 


Optical_Depth_Ratio_Small_Ocean_055micron 

(X,Y,2): best, average 


Optical_Depth_by .models 

(X,Y,9): 9 models 


Number_Pixels_Used_Ocean 

(X,Y) 


Mean_Reflectance_Ocean 

(X,Y,7); 0.47, 0.55, 0.66, 
1.24, 1.63, 2.13 |am 

CD 

00 

o 

STD_Reflectance_Ocean 

(X,Y,7); 0.47, 0.55, 0.66, 
1.24, 1.63, 2.13 [am 

co" 

00 

o 

Quality _Assurance_Ocean 

(X,Y); Packed byte 



X, Y refers to a 2-dimensional array along and across the swath. Some parameters have a third dimension. 
A dimension of “7” refers to the 7 retrieved ocean wavelengths, listed. A dimension of “2” refers to either the 
solution with minimum fitting error (best) or the average of all solutions with fitting error less than 3% 
(average) (Remer et al., 2005). A dimension of “9” refers to the 9 models, 4 fine mode and 5 coarse mode, 
used in the retrieval (Remer et al., 2005). 
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Table 3. Land parameters of the MOD/MYD04_3K file and the parameter’s dimension. 


Parameter 

dimension 

Aerosol_Type_Land 

(X,Y) 

Fitting Error_Land 

(X,Y) 

Su rface_Ref lectance_Land 

(X,Y,3); 0.47, 0.66, 2.13 pm 

Corrected_Optical_Depth_Land 

(X,Y,3*); 0.47, 0.55, 0.66 pm 

Corrected_Optical_Depth_Landwav2p1 

(X,Y) 

Optical_Depth_Ratio_Small_Land 

(X,Y) 

Number_Pixels_Used_Land 

(X,Y) 

Mean_Reflectance_Land 

(X,Y,7); 0.46, 0.55, 0.66, 0.86, 1.24, 1.63, 2.13 pm 

STD_Reflectance_Land 

(X,Y,7); 0.46, 0.55, 0.66, 0.86, 1.24, 1.63, 2.13 pm 

Mass_Concentration_Land 

(X,Y) 

Aerosol_Cloud_Fraction_Land 

(X,Y) 

Quality _Assurance_Land 

(X,Y) 

Topographic_Altitude_Land 

(X,Y) 


X,Y refers to a 2-dimensional array along and across the swath. 

Some parameters have a third dimension. 

A dimension of “3” refers to the 3 wavelengths of input reflectances used in the retrieval, listed. 
A dimension of “3*” refers to the 3 retrieved wavelengths over land, listed. 

A dimension of “7” refers to the 7 solar wavelengths, listed. 
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Table 4. Parameters from MODIS-AERONET collocation validation: means of each data set, 
correlation coefficient, regression slope and offset, number of collocations, percent within ex- 
pected error*, percent above expected error* and percent below expected error*, for 10 km abd 
3 km products, separated into land and ocean retrievals. The last rows show MODIS-MAN col- 
locations over ocean. 



AERONET 
or MAN Mean 

MODIS 

Mean 

R 

Slope 

offset 

N 

% 

EE 

% 

above EE 

% 

below EE 

land 10 km 

0.158 

0.172 

0.87 

0.93 

0.02 

3298 

71 

18 

11 

land 3 km 

0.153 

0.200 

0.84 

1.04 

0.04 

3283 

62 

31 

7 

ocean 1 0 km 

0.143 

0.148 

0.94 

0.95 

0.01 

1100 

69 

20 

11 

ocean 3 km 

0.139 

0.156 

0.90 

0.96 

0.02 

697 

69 

24 

7 

Ocean 1 0 km MAN 

0.072 

(median = 0.054) 

0.094 

(median = 0.064) 

0.93 

1.04 

0.02 

37 

65 

35 

0 

Ocean 3 km MAN 

0.077 

(median = 0.051) 

0.101 

(median = 0.059) 

0.94 

1.21 

0.01 

37 

68 

30 

2 


* Expected error for land is ±0.05 ±0.15 AOD, and for ocean it is ±0.03 ± 0.05 AOD. 
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} 



— 3 km— ► 

6x6 pixels 


* 10 km ■ » 

20 x 20 pixels 


Fig- 1- Illustration of the organization of the MODIS pixels into retrieval boxes for (left) the 1 0 km 
product consisting of 20 x 20 0.5 km pixels within the magenta square and (right) the 3 km 
product consisting of 6 x 6 0.5 pixels within the red square. The small blue squares represent 
the 0.5 km pixels. The white rectangles represent pixels identified as cloudy. The 3 km retrieval 
box is independent of the 10 km box, and is not a subset. Here it is shown enlarged. 
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Fig. 2. Flowchart illustrating the different paths of the 10 km (red) and 3 km (blue) retrievals. 
The procedures appearing in the black outlined boxes are common to both algorithms. 
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AOD at 550 nm 

- 0.05 0.05 0.15 0.25 0.35 0.45 0.55 0.65 0.75 


Fig. 3. Aerosol optical depth at 550 nm retrieved from the 15 July 2010 Aqua-MODIS radi- 
ances using the Collection 6 MODIS Dark Target aerosol algorithm. Left: the product at 10 km 
resolution. Right: the product at 3 km resolution. The situation is a moderate dust event over 
the Mediterranean Sea off the costs of Tunisia and Libya. The 3 km retrieval produces values 
closer to the coastline and to the islands. 
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Fig. 4. Same as Fig. 3 but for a situation over Canada where the 3 km product better resolves 
the plume from active wildfires. Note the additional red pixels between the clouds in the upper 
left corner of the image for the 3 km panel. 
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AOD at 550 nm 


•0.05 0.05 0.15 0.25 0.35 0.45 0.55 0.65 0.75 

Fig. 5. Same as Fig. 3, but for 12 January 2010 during a pollution episode in China. Here 
the 3 km algorithm is able to make retrievals over a much broader region. The AOD interpo- 
lated to 0.55 pm from the only AERONET station in the image (Hong_Kong_Polyll) is 0.38. The 
3 km retrieval there is 0.45, while there is no 10 km retrieval available at that spot during this 
overpass. 
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Fig. 6. Same as in Figure 3 but for a situation over California where the 3 km product introduces 
artificial noise over an urban area that the 10 km product avoids. 
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Fig. 7. Difference between 3 km daily global mean aerosol optical depth (AOD) at 550 nm and 
the same at 10 km for over ocean (left) and over land (right). Positive differences indicate that 
the 3 km AOD is higher than the 10 km AOD. 
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Fig. 8. Daily global mean aerosol optical depth (AOD) at 550 nm for the 93 days of the Januarys, 
and 93 days of the Julys 2003, 2008 and 201 0 calculated from the 3 km product plotted against 
the same calculated using the 10 km product. The left plot is for ocean, and the right plot is for 
land. 
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Fig. 9. Frequency histograms of 10 km and 3 km AOD retrievals, constructed from global re- 
trievals from three Januarys (2003, 2008 and 2010) in the left column and from three Julys 
(same years) in the right column. Land in the top row and ocean in the bottom. 
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Fig. 10. Bin by bin differences between the frequency histograms of Fig. 9 (3km-10km). 
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Fig. 11. Binned scatter plot of MODIS AOD at 550 nm retrieved from over land against 
AERONET observed AOD. MODIS 10 km product on the left and 3 km product on the right. 
Each bin represents 50 collocations. AERONET AOD has been interpolated to this wavelength 
to match MODIS. Error bars represent ±1 standard deviation in the bin. The green line is the 
1:1 line. The red line is the linear regression through the full database of 3298 and 3283 col- 
locations for 10 km and 3 km, respectively. The dashed lines represent expected error of the 
MODIS 10 km retrieval. Regression statistics are given in Table 4. 
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Fig. 1 2. Same as Fig. 1 1 , but for over ocean retrieval. Here, the linear regression line is in light 
blue and the 1 :1 line is in black. 
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Fig. 13. Binned differences of MODIS - AERONET AOD 550 nm plotted as a function of 
AERONET AOD. Land is on the left and ocean on the right. 3km products are denoted by 
dashed lines and 10 km products by solid lines. The x-axis has been truncated at AOD = 0.6 to 
emphasize the lower AOD bins. 
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MAN AOD 550 nm MAN AOD 550 nm 


Fig. 14. Scatter plot of MODIS 10 km and 3 km retrievals of AOD at 550 nm plotted against 
collocated data from the Maritime Aerosol Network (MAN). 37 collocations were identified at 
each resolution in the 6 months undergoing analysis, although they are not necessarily the 
same 37 days in each plot. 
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